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ABSTRACT 

We construct test-particle orbits and simple N-body models that match the properties of the 
giant stellar stream observed to the south of M31, using the model of MSl's potential de- 
rived in the companion paper by Geehan et al. (2005). We introduce a simple approximation 
to account for the difference in position between the stream and the orbit of the progenitor; 
this significantly affects the best-fitting orbits. The progenitor orbits we derive have orbital 
apocenter ^ 60 kpc and pericenter ^ 3 kpc, though these quantities vary somewhat with the 
current orbital phase of the progenitor which is as yet unknown. Our best combined fit to the 
stream and galaxy properties implies a mass within 125 kpc of M31 of (7.4± 1.2) x 10^^ Mq. 
Based on its length, width, luminosity, and velocity dispersion, we conclude that the stream 
originates from a progenitor satellite with mass Mg ^ 10® ™d most modest amounts 
of dark matter; the estimate of Mg is again correlated with the phase of the progenitor M3 1 
displays a large number of faint features in its inner halo which may be progenitors or contin- 
uations of the stream. While the orbital fits are not constrained enough for us to conclusively 
identify the progenitor, we can identify several plausible candidates, of which a feature in the 
planetary nebula distribution found by Merrett et al. is the most plausible, and rule out several 
others. We make predictions for the kinematic properties of the successful candidates. These 
may aid in observational identification of the progenitor object, which would greatly constrain 
the allowed models of the stream. 

Key words: galaxies: individual: M31 - galaxies: interactions - galaxies: kinematics and 
dynamics 



1 INTRODUCTION the extended tidal tails of the Sagittarius dw arf galaxy that wrap 

around the Milky Way (Totten & Irwin 1998), though its azimuthal 
The eiant southern stream observed near the Andromeda galaxy . ^ j »ni ■ i ha ■ j-cr ■ . 

" o J extent around M31 is much smaller. A major diiierence is that we 

(M31) provides us with an extraordinary opportunity to study the 



disruption by a large disk galaxy of one of its satellite galaxies, 
and thereby obtain a sensitive probe of the dynamics in M3 1 . The 
stream was originally discovered through counts of red giant stars 



do not know if the progenitor of the Andromeda stream has sur- 
vived and if so where it is, although several faint features in the 
inner halo of M31 have been suggested as possible progenitors. 



in the M31 halo (Ibata et al. 2001a). Later extensions of this survey Since the three-dimensional path of the stream is now mea- 

( Ferguson et al. 2002; McCo nnachie et alj 120 03 ; Ferguson et al. sured along with one velocity component, we have nearly com- 



12004) outlined the extent of the stream from a projected distance plete information on the path of the system in phase space. If we 

of 60 kpc from the center of M3 1 on the southern side, to a pos- combine the stream data with other observations of M3 1 such as 

sible extension o n the northern side of M31. Fits to the tip of the its rotation curve, we should be able to obtain a direct estimate 

red giant branch iMcConnachie et al. 2003) provided line-of-sight of the gravitational potential in the M31 halo, as well as follow 

distances as a function of position along the stre am. Finally, spec- the path of the stream beyond the visible southern po rtion. Cur- 

troscopy of stars along t he path of the stream ilbata et al]|2004t rent w ork includes some first steps in this direction. iFont et all 

iGuhathakurta et alj2005h has resolved the stream as a distinct kine- i2005h computed orbits in a bulge-disk-halo model. These authors 

matic component as well, and has also provided absorption-line es- found rough consistency of their potential with the stream, but 

timates of the metallicity of the stream stars. The stream resembles they did not optimize the agreement in detail, focusing instead 

on the physica l implications of t he width and velocity dispersion 
of the stream. Ilbata et all i2004i) tested several potentials, deriv- 

* E-mail: fardal@fcraol.astro.umass.edu, babul@uvic.ca, jgee- ing a weakly model-dependent limit on the mass within 125 kpc. 

han@uvic.ca, raja@ucolick.org They supplied model parameters only for their simple spherical po- 
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tentials, which have the disadvantage that they do not match the 
M31 rotation curve. They used a more sophisticated bulge-disk- 
halo model as well, but this is not fully specified in the paper. 

An issue that deserves closer attention is the slight difference 
between th e locations of the orbit and the stream. Both I bata et al. 
J2004 and lFontetal]j2005h assume for simplicity that the stream 
follows the orbit. However, the stream exists precisely because its 
constituent stars are on slightly different orbits. The systematic vari- 
ation of orbital energy along the stream, which is accounted for 
implicitly by N-body simulations but not by simple orbital integra- 
tions, may affect the estimates of the M31 potential and the trajec- 
tory of the progenitor. Hence it is worth building upon the earlier 
work to construct a more specific and more accurate model of the 
stream. 

In lGeehan et alj I200I (Paper I), we presented a simple ana- 
lytic bulge-disk-halo model for the potential of M31, fitted to ob- 
served tracers of the potential ranging from bulge stars to satel- 
lites in the outer halo, and we quantified the error on the model 
parameters. In this paper, we use this fit to derive simple models 
of the southern stream, and discuss their physical implications. We 
start off with a subset of the observational data and then gradually 
add in more constraints. In § 2, we consider the orbital informa- 
tion alone. We integrate orbits in our M31 potential to obtain best 
values and confidence intervals for the orbital parameters, and fur- 
ther refine the galactic parameters compared to Paper I. We intro- 
duce an approximate method to calculate the phase space location 
of the stream given the orbit of the progenitor, test it with N-body 
simulations, and incorporate it within the fitting procedure; as a re- 
sult the allowed orbits are changed considerably. In § 3, we add 
in the constraints that the observed length, width, velocity disper- 
sion, and luminosity of the stream put on the properties and current 
location of the progenitor satellite, using N-body simulations that 
follow our best-fitting orbit. § 4 examines the forward continuation 
of the stream, including the progenitor location, and its relation to 
the orbital characteristics of the stream. Among the numerous faint 
features in the inner halo of M31, we can identify some that are 
plausibly related to the stream. Finally, § 5 presents our conclu- 
sions. 



2 ORBITAL FITS 

2.1 Summary of the observations 

Though it was discovered only recently, the variety of observa- 
tional information about the southern stream already approaches 
that which is available for the Sagittarius stream in our own galaxy. 
We review here what is currently known about the configuration of 
M3 1 and the southern stream. 

We look at the M3 1 galaxy from below, as is evident from im- 
ages of the disk's dust lanes projected on the bulge, and the galaxy 
spins counterclockwise on the sky. We take the inclination of M3 1 
to be 77°, and the position angle to be 37°. We use a distance to 
M31of2g = 784±24kpc(Stanek&Garnavich 1998), and a mean 
radial velocity of Vgz — —300 ± 4kms~^ ide Vaucouleurs et alJ 
I1991I) . 

Because the M31-stream system covers such a large region 
on the sky, its transverse velocity is one of the primary uncertain- 
ties in its internal kinematics, which is a situation rarely encoun- 
tered in astronomy. Since M31 and the Milky Way are the pri- 
mary mass concentrations in the Local Group, it is traditionally 
assumed that they are falling in toward each other for the first 
time, so that the transverse velocity is much smaller than the ra- 



dial velocity. We prefer to use the results of IPeebles et alj i200lh . 
who have modeled the Local Group dynamics in detail. They 
found solutions for M3rs transverse velocity falling mostly into 
two clumps in opposing directions, with velocities in supergalac- 
tic coordinates I and b of either {vgi,Vgb) = (—40, 120) km s^^ 
or (30, — 140)kms~^. The eastward and northward velocities are 
then {vg^,Vgy) = (-130, 20)kms"^ or (150,0)kms"\ The 
dominance of the east-west component stems from the direction of 
the local filament, along which M3 1 tends to move in the Peebles 
et al. model. 

The giant southern stream extends to the SE of the galaxy for 
approximately 5°. The total extent of the stream is not well known; 
at its southern end it nuiy be visibly petering out close to Field 1 of 
McConnachie et al. (2003), as seen in the larger survey presented 
in Fer guson et al. 1,2004.) . and continuing northwards it is eventu- 
ally lost to confusion against the M31 d isk. The stream luminosity 
is ~ 3.4 X 10^ Lq ^Ibata et aljEoOI j) . suggesting a stellar mass 
~ 2.4x 10* Mq if we assume M/Ly ~ 7. The mean me tallicity of 
the stream is [Fe/H] ~ —0.5 iGuhathakurta et ai]l2005h . suggest- 
ing the progenito r had a much larger total luminosity of ~ lO'' Lq 
( iFont et alj |2QQA. McConnachie et al. suggest the stream contin- 
ues at least 1.5° to the NW on the near side of M31, but this is 
questionable for several reasons: t he purported extensio n is not a 
distinct surface brightness feature jFerguson et aljEoo4) . there is 
as yet no spectroscopic confir mation of a stream component, and 
as shown in Ibata et al] J2004 it is difficult to derive an orbit that 
fits this feature. We therefore consider this "northern extension" to 
be just one of many faint features that could possibly be extensions 
of the stream. In contrast, the reality of the southern portion is in no 
doubt, due to the strong enhancement in star counts and observa- 
tion of a distinct kinematic component in five spectroscopic fields 
along the stream. 

Field coordinates for the published observations of the stream 
are listed in Table Q along with distances and velocities of the 
stream. The positions show the stream is a roughly linear structure, 
at least in projection. The field "a3" of Guhathakurta et al. 1 200^ 
overlap s the stream only by chan ce, and is not centered on the 
stream. iMcConnachie et alJ^2003^ chose their fields to lie approx- 
imately along the stream center. However, some of these positions 
were chosen by extrapolation from the inner part of the stream. 
In fact, their fig3, showing a histogram of stars as a function of 
transverse position, suggests the central ridge of the stream is off- 
set fairly uniformly by ~ 0.15° from the field centers, towards the 
SE. We use the slope of Fields 1-8 to define rotated sky coordinates 
m and n, where m = 0.504^ — 0.864r; increases along the stream 
(roughly toward the SE) and n = -0.864^ — 0.5047; increases 
across the stream (roughly toward the SW). We then shift the mean 
transverse position of the stream relative to the field centers by 
0.15°. This shift is small compared to the width of the stream, ap- 
proximately 0.5° wide full width at half-maximum (FWHM) on 
the sky or about 7 kpc in distance, as seen in the same figure. We 
assign an uncertainty of 0.15° to the stream's angular position. 

D istances to the M31 stream fields are given by 
IMcCo nnachie et al. (2003) and reproduced in Table Q These 
show the stream lies further away than M31, and the stream 
follows a roughly radial trajectory. We exclude the fields on the 
northern side, and exclude the field nearest to M3I (Field 8) as well 
because of possible contamination by disk stars. McConnachie 
et al. estimate a total error of 20 kpc on their distances, but this 
is composed of ~ 12 kpc random error and ~ 16 kpc systematic 
error from the photometric zero-point and the stellar population 
uncertainties. Orbits that have a steeper or shallower gradient 
of distance along the stream than the observed gradient, as in 
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Table 2. Parameter choices for some orbits discussed in the text. The first six parameters give the initial position and velocity of the test particle (the orbit i s 
calculated both forward and backward of this point). Az is the fitted systematic offset between our distance calibration and that off McConnachie et alJi20o3) . 
The halo parameter ff^ determines the strength of the galaxy potential. The final parameter gives the value of reduced from our fit. Initial conditions le, 
2e, 3e, and 4e give the same orbits as 1,2, 3, and 4, except they begin at an earlier point on the orbit. 



Orbit 


x{0) 


2/(0) 


z{0) 


v^{0) 


Vy{0) 




Az 


fh 


Fp 




comments 






kpc 


kpc 


kpc 


km s~^ 


km s~^ 


km s^^ 


kpc 












1 


5.19 


—20 


38.0 


—61.2 


132.9 


—214.4 


—4.0 


1.14 




0.16 


fit to orbit, variable halo 






il.l2 




±8.8 


±9.4 


±21.9 


±6.7 


±12.8 




±0.35 








le 


—52.7 


42.6 


-28.2 


126.5 


—94.2 


53.9 


—4.0 


1.14 




0.16 


same orbit, earlier stall 




2 


5.60 


—20 


18.4 


—67.9 


145.2 


—209.5 


—21.42 


1.14 


1.00 


0.51 


fit to stream, vai' halo, Fp 


= 1.0 




±1.02 




±2.7 


±4.0 


±13.7 


±7.3 


±6.4 


±0.37 










2e 


—8.8 


— 19.1 


—56.4 


4.9 


44.8 


42.4 


—21.42 


1.14 


1.00 


0.51 


same orbit, earlier 




3 


7.04 


—20 


23.6 


—51.9 


85.2 


— 155.1 


— 18.70 


1.16 


1.50 


0.48 


fit to stream, var halo, Fp 


= 1.5 




±0.66 




±5.8 


±6.5 


±14.3 


±8.5 


±11.4 


±0.40 










3e 


-5.5 


-11.8 


-42.1 


5.9 


45.7 


53.2 


-18.70 


1.16 


1.50 


0.48 


same orbit, earlier 




4 


5.15 


-15 


16.7 


-59.4 


102.3 


-175.0 


-19.88 


1.10 


2.00 


0.46 


fit to stream, vai' halo, Fp 


= 2.0 




±0.51 




±2.8 


±4.1 


±12.0 


±7.8 


±8.0 


±0.38 


±0.00 








4e 


-2.2 


-7.5 


-32.2 


4.4 


52.8 


81.1 


-19.88 


1.10 


2.00 


0.46 


same orbit, earlier 




5 


8.54 


-20 


29.4 


-7.2 


6.4 


-57.7 


-15.55 


1.02 


2.02 


0.97 


fit to M32 






±0.60 




±2.4 


±17.9 


±28.1 


±54.5 


±5.1 


±0.32 


±0.01 








6 


5.78 


-20 


25.7 


-62.4 


118.9 


-196.7 


-12.86 


1.47 


1.12 


0.57 


fit to Northern Spur 






±0.93 




±3.3 


±7.4 


±16.8 


±9.3 


±8.1 


±0.46 


±0.03 








7 


6.96 


-20 


25.1 


-51.0 


83.2 


-154.5 


-16.20 


1.18 


1.50 


0.47 


fit to And NE 






±0.71 




±2.8 


±9.4 


±17.9 


±12.7 


±6.4 


±0.40 


±0.13 








8 


5.25 


-15 


24.2 


-63.6 


136.4 


-220.9 


-7.07 


2.21 


1.29 


1.00 


fit to Merrett NE PNe 






±0.52 




±3.2 


±8.0 


±19.2 


±7.7 


±9.1 


±0.68 


±0.04 









Table 1. Kinematic data for the southern stream, in units where M3 1 is at 
the center A dash indicates that there is no data for that field. The angular 
positions g and t) are those of the field centers; field "a3" is certainly offset 
from the position of the stream, while fields 1-8 may be as well. 



Field 


? (deg) 


V (deg) 


d' (kpc) 


v'j. (km s ^) 


a3 


+1.077 


-2.021 


66 


-158 


1 


+2.015 


-3.965 


102 





2 


+1.745 


-3.525 


93 


-50 


3 


+1.483 


-3.087 


76 




4 


+1.226 


-2.653 


71 




5 


+0.969 


-2.264 


56 




6 


+0.717 


-1.768 


52 


-180 


7 


+0.467 


-1.327 


45 




8 


+0.219 


-0.886 


-4 


-300 


12 


-0.731 


+0.891 


-45 




13 


-0.963 


+1.342 


-26 





fig. 4c of'ibata et alj j2004 or fig. 2a of lFontetalJ llOOdt) . can 
thus be ruled out even if they lie within the total error bars of 
McConnachie et al. To treat this in our orbital fitting, we take 
the observed distances to be the McConnachie values plus a free 
parameter, Az. Assuming the uncertainty on this systematic shift is 
aAz ~ 16 kpc, we then include a term (zXz/ctaz) in the total x ■ 
We set the random error on each distance point equal to 12 kpc. 

The mean velocity at various points al ong the stream has been 
measured in four fields along the stream by Ibata et al. 72004"), and 
in the single stream field "a3" of Guhathakurta et al. (2005). This 
reveals that the stream is falling in toward M31 from behind. We 
again omit the field closest to M3 1 (Field 8 of Ibata et al.) because 
of contamination concerns. iGuhathakurta et al. 1 200^ e stimate the 
random error on their mean velocity to be 4 km s~^ . Ilbata et'al] 



yOOJ) do not quote errors on their mean velocities, but they are 
probably comparable since they estimate the velocity dispersion it- 
self to be only llkms"^. However, the true uncertainty includes 
systematic effects such as the 4kms^^ error in MBl's radial ve- 
locity, and contamination by the possibly lumpy velocity distribu- 
tion of the general halo population. Because of these concerns, we 
(somewhat arbitrarily) raise the estimate of the velocity error to 
8kms"^ 

Our orbits are computed in a rectangular right-handed coordi- 
nate system {x, y, z), centered on M3 1 and oriented with respect to 
the sky so that x points east, y points north, and z points away from 
us. However, the observed quantities are not obtained in a rectangu- 
lar coordinate system, because M3 1 is so close that the system has 
a significant width and depth on the sky. Insufficient information is 
available to correct all the observations to rectangular coordinates, 
so we convert the orbits to observed coordinates instead. The posi- 
tional coordinates ^ and rj are standard coordinates on the tangent 
plane. We obtain the east-west coordinate ^, the north-south coordi- 
nate r], the distance d, and the radial velocity Vr using the following 
formulae: 

e = x/iz, + z) (1) 

V = y/{zg + z) (2) 
d = d' + Zg ^ [x'^ + y"^ + {zg + z)'^Y''^ (3) 

Vr = [x{Vj:+Vgj:) + y{Vy+Vgy) + {Zg + z){Vz+Vgz)](4-) 

= V'^ + Vgz 

In our plots we use the relative distance d' (which we call the 
"depth") and the relative velocity v'^., found by subtracting the mean 
distance and radial velocity of M31. These equations could also 
be expanded in orders of M31-centric quantities, since in practice 
only the first-order correction terms are significant. For the two 
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solutions for the transverse velocity of M31 {vgx,Vgy) discussed 
above, the effect of MSl's transverse velocity on v,. is quite mod- 
est, only about ±5kms^^ for Field 1 which is the most strongly 
affected field. This is partly because these velocity solutions point 
almost perpendicular to the stream. In contrast, a maxi mum veloc- 
ity of ±300 km s~^ along the stream, as allowed for bv llbata et alJ 
yOOJ), would have added a significant uncertainty of ~ 30 km s~^. 
We set Vgx = Vgy = from now on. 



2.2 Fits to stream directly 

We now try to find orbits that fit the stream observations. At 
first, for simplicity, w e assume th at the stre am follows the or - 
bit, as in the work of llbata et J1 ^2004) and iFont et'ZI j2005h . 
In Paper 1, we demonstrated that even with this assumption 
our orbits can still differ significantly from those in the cited 
papers, because of the different radial mass distribution and 
the flattening of the disk potential. The possible flattening of 
the halo component is a contentious issue e ven in the case of 
the Sagittarius stream in ou r own Galaxy jibataetalJ I2001bl: 
Marti'nez-Delgado et al. '2004'; 'John ston. Law. & Maiewskil l20~ ; 
Law, Johnston, & Majewski 2005), which is more suitable for de- 
tecting the effects of a non-spherical potential due to its much larger 
azimuthal extent. In any case the disk contributes significantly to 
the overall flattening in the Milky Way, and thus probably in M3 1 
as well. In this paper we will simply assume the halo is spherical, 
but this is certainly an issue to examine in the future. 

We ignore the effect of dynamical friction on the satellite orbit. 
This is partly for convenience, as the strength of dynamical friction 
depends on the unknown mass and size of the progenitor, and thus 
complicates the estimates we make of these quantities in § 3 below. 
More importantly, we expect the effect of dynamical friction on the 
energy to be fairly small (< 10% per orbit) unless the progenitor 
satellite mass Ms is quite large (Ms ^ 10^ Mq). Our simula- 
tions typically have Ms ~ 10^ Mq at the start of the run, and the 
effects of tidal stripping and gravitational disk shocking further re- 
duce this mass during the run. Dynamical friction due to the disk 
itself is lik ely important only for satellites with orbits in the plane 
of the disk iTavlor & Babul'200l'). Furthermore, for a highly radial 
orbit such as that required by the southern stream, the energy loss 
is almost entirely due to the passage through pericenter. The energy 
within the stream, which is observed over less than one radial pe- 
riod, is thus affected very little by the neglect of dynamical friction. 
In § 4, we speculate on the location of the progenitor allowing for 
the possibility that it has passed through pericenter once again. The 
neglect of dynamical friction there may induce a small error, the 
size of which depends on how much of the mass has managed to 
remain as an intact object during this second passage. Accounting 
for this would greatly complicate the analysis, so we leave treat- 
ment of dynamical friction to future work. 

We also assume a static potential, ignoring any interaction 
with M3rs satellite galaxies or recoil and/or distortion of M3rs 
potential. The latter assumption is justified because of the large 
mass of M31 (~ 10^^ Mq in the bulge and disk alone), compared 
to the typical estimated mass of the progenitor (~ 10^ Mq) that we 
estimate below. With these assumptions, the zero-point of the orbit 
is arbitrary. We absorb this freedom by setting the orbit to begin in 
a fixed plane, usually j/(0) — —20 kpc. 

Paper I showed that there is considerable uncertainty in the 
potential of M3 1 , especially in the halo regions where the stream 
lies. This uncertainty is mostly due to two degeneracies between 
the model parameters, namely a tradeoff between the disk and halo 



components and between the scale radius and density parameter 
5c of the halo. Much of the formally allowed region, however, im- 
plies unlikely values for the physical parameters like the halo con- 
centration, baryonic fraction, and disk mass-to-light ratio; thus the 
likely region of parameter space is much smaller. Because the de- 
tailed properties of the disk and bulge are not probed by the stream, 
we limit our exploration of parameter space to a one-dimensional 
path that varies the strength of the force at large radii. We set the 
disk surface density to a fixed value Eo = 4 x 10* Mq kpc~^; 
the bulge and disk scale lengths rb = 0.61 kpc and rd = 5.4 kpc 
are already fixed by our fit to the surface brightness. (These param- 
eters are defined completely in Paper I.) We then find the best fit 
values of the bulge mass Mi, and the halo density parameter 5c as a 
function of the halo scale radius r^, storing them in a lookup table 
along with the value of x^- In other words, we are probing a verti- 
cal path through the contours in fig. 6 of Paper I. For convenience, 
we describe the location along this path with a "halo parameter" 
fh = rh/r^^\ where r^*^' = 9.0 kpc is the best fit scale radius for 
the chosen Eo in our flattened disk model. To within a few percent, 
the best-fitting bulge mass and halo density parameter along this 
path are given by 

Mb = (2.71 + 0.435it - Q.lQlu + 0.0229u^) x 10^° Mq , 
5c = exp(12.46 - 1.928u + 0.1434u^) , 

with u = In fh - The mass within 125 kpc, approximately the max- 
imum radius of the stream, scales roughly as M125 = 6.68 x 
10"/^-^^ Mq. The allowed deviation of Sc from its best-fitting 
value is very small once the other parameters are specified, due to 
the small error bars on the rotation curve, so we do not try to adjust 
it from its best fit value to vary the halo potential further. 

This gives seven free parameters to consider: x(0), z(0), 
Vx{Q), Vy{Q), v^{Q), Az, and fh. x{0) and the ratio Vx{0)/vy{0) 
are fairly tightly constrained by the positional data, while Vz{0) 
is constrained by the spectroscopic data. This leaves z{0) and 
Vx (0) /vz (0) as the most uncertain parameters, since they involve 
the distance measurements. 

We use a simple Runge-Kutta integration procedure to trace 
the orbit, since extremely high accuracy over many orbital periods 
is not required. We use a Levenberg-Marquardt routine to opti- 
mize the fit to the stream. 

Our best-fitting orbit (labeled orbit 1 in Table |2j is shown as 
the solid line in Fig. Q We also show the same orbit calculated 
with a constant conversion of transverse distance to angle, using the 
distance to M31 of 784 kpc. When this is done, the stream moves 
outward significantly in ^ and rj because it lies further away than 
M31. v'^ also changes by a small amount ( < 5 km s^^) due to the 
variation in the radial vector on the sky. 

Our best-fitting orbit fits the data quite well; in fact, the ~ 
2.2 for 14 degrees of freedom is almost too good. (In part, this is 
because the positional data points are derived from the field cen- 
ters rather than calibrated individually, and thus have effectively 
been smoothed.) The best-fitting parameter values and their for- 
mal errors from the covariance matrix are given in Table |2| The 
radial period is 1.7 Gyr, apocenter is at 119 kpc, and pericen- 
ter is at 2.2 kpc. The orbit reaches apocenter near the most dis- 
tant field (Field 1); this happens because the velocity relative to 
M31 is close to zero, requiring turnaround there. If we include 
only the observational constraints from the stream, the preferred 
halo parameter is fh = 1.8 ± 1.1, corresponding to a mass of 
i\fi25 — (9.1 ± 3.1) X 10^^ ^0- This suggests a more massive 
halo than in our preferred galaxy model, which has fh ~ 1 .0 ± 0.4, 
or M125 = (6.7 ± 1.3) x 10^^ Mq. The observational constraints 
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Figure 1. Best-fitting orbit assuming the stream follows the orbit, using the flattened disk potential. The orbit is shown with purple solid lines. The orbit 
calculated using a rectilinear coordinate system is shown with red dotted lines (see discussion in the text). Observational results are shown with black symbols. 
(Upper left): North-south coordinate t) versus east-west coordinate ^. The direction vectors of the stream coordinates m and n are shown in the plot. Crosses 
indicate the stream positions derived from the field coordinates of McConnachie et al. i2003), offset slightly as discussed in the text; the fields are labeled 
according to the numbering in that paper. (Upper right): Transverse stream coordinate n, versus distance along stream m, with the vertical scale expanded for 
clarity. Diamonds indicate the same stream positions as in the pr evious panel, again labeled by field number (lower left): Depth d' relative to M31 versus m. 
Diamonds indicate the observational constraints of McConnachi e et alj i2003h . as adjusted by the systematic shift Az from our fit; open circles give the depth 
values with the original calibration. The leftmost diamond is left empty as a reminder we do not use this point in the fit, due to concerns over contamination 
by the disk. The error bar in the bottom right comer indica tes t he uncertainty on Az. (Low er right): Radial velocity ri^ relative to M31, versus m. Diamonds 
indicate the observational consfi'aints of llbata et alj l2004l) and lGuhathakurta et alj i2005ft . The leftmost symbol is again left empty as a reminder we do not 
use this point in the fit. 



on the galaxy fit are stronger than those on the stream, however, so 
when we combine the stream data with the galaxy fit the preferred 
halo scaling parameter is raised only slightly: fh — 1.1 ± 0.3, 
or Mi25 = (7.2 ± 1.0) X lO" Mq. The fits with and without 
the galaxy term are both consistent with Iba ta et aL who 
found a Icr confidence interval of M125 = (6.2-10.0) x 10^^ Mq. 
The orbit is less strongly radial than in Ibata et al., which has 
Rapo ~ 1 25 kpc and R 



1.8 kpc. It is more radial than 



the orbits in lFont et alj i200^ . which have Rapo ~ 110 kpc and 



Rp 



4 kpc. The direction of their second loop is also signif- 



icantly different. The differences between these various orbits are 
induced both by different assumptions about the potential of M31 
and by differences in the chosen fitting technique and treatment of 
the observational data. 



2.3 Stream-orbit tilt 

The next step is to incorporate a phase space tilt between the 
orbit and the stream. When a satellite is tidally disrupted, stars 
on different orbits arrange themselves primarily as a function of 
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their orbital energy onto a stream that nearly follows the orbit 
jjohnston. Sackett. & Bulloc3l200lh . The "nearly" is important, 
though: the stream is tilted relative to the progenitor's orbit, both 
in 3-d real space and 6-d phase space. For an object like the 
Andromeda stream which is on a highly radial orbit, the stream 
stars are lost from the progenitor almost entirely at pericenter. If 
the progenitor is disrupted over the course of several radial pe- 
riods, a stream results from each pericentric passage, with the 
youngest stream hav ing the largest tilt relative to the orbit. (Fig. 1 in 
iLaw. Johnston. &Maiewski.2005, displays these effects for a model 
of the Sagittarius stream.) The length of the stream increases both 
with time and with the scatter in the energy; this scatter depends in 
turn on the satellite mass. The phase space tilt between the orbit and 
the stream, however, decreases with time and to first order is inde- 
pendent of the satellite mass. Hence, to predict the stream location 
given the progenitor's orbit, we need only to know when the stream 
stars were removed from the progenitor, and not the mass or other 
properties of the progenitor. In theory this could include pericen- 
tric passages before the last one, but this is unlikely beca use of the 
narro wness and high surface brightness of the stream iFont et alj 
l2005h . (Note: when we refer to the "last" pericentric passage, we 
mean the last one experienced by the southern stream stars, unless 
otherwise specified. The progenitor itself may already have passed 
through pericenter again.) 

The current location of the progenitor along the orbit thus be- 
comes an important variable. This is currently unknown, but we 
can at least put some constraints on it. The progenitor is unlikely to 
trail the southern stream in its orbit, since there would be another 
stream behind it which is not observed. The progenitor is not ob- 
served along the stream, although it could be located there if it is 
almost entirely disrupted. However, the lumi nosity of the proge n- 
itor inferred from the stream's metallicity bv lFont et alj i2005h is 
much larger than the total luminosity of the stream, suggesting this 
is rather unlikely. If the progenitor lies too far ahead of the stream 
in its orbit, on the other hand, the corresponding stream length and 
hence the implied scatter in the orbital energies may require a im- 
probably large progenitor mass, as discussed below. 

Suggested progenitors for the stream include the And VIII 
concentration, the Northern Spur, and a high-metallicity feature on 
the east side of the disk. All lie near the M3 1 disk, and thus would 
have to precede the stream in its orbit. If one of these was the source 
of the stream, one might expect that if it had passed the next peri- 
center and was heading inbound again, there would be an obvious 
continuation of the stream pointing outwards. The lack of such a 
feature is a weak argument that the orbital phase of the progeni- 
tor is < 1.5. It is also possible that the progenitor lies somewhat 
ahead of the southern stream in the vicinity of pericenter, but is 
projected against M3I and has escaped detection completely. These 
arguments suggest the stream stars were most likely unbound from 
the satellite in the range 0.8-1.5 radial periods ago, at the last peri- 
center passage. 

To fit the orbit, many orbital integrations are required, so com- 
puting the tilt between the stream and the orbit using an N-body 
simulation on each integration would be extremely cumbersome. 
We need a simple method to infer the position and velocity of the 
stream, given the orbit of its prog enitor. Our model builds on the 
ideas presented in IjohnstorjITgga) . That paper shows the azimuthal 
phase at apocenter does not vary strongly with orbital energy, i.e., 
along the stream. From this observation, it infers that the azimuthal 
phases of the different stream orbits are self-similar. That is, the az- 
imuthal phase ^ofa given orbit in the stream as a function of time 
t past pericenter is simply ^{t) = *o + '^cen{t/T + to), where 
Vl/cen is the phase of the central orbit, the time scaling factor r de- 



pends on the orbital energy, and to and ^Po set the time and azimuth 
at the point where the progenitor is disrupted, creating the stream. 

We extend these ideas by suggesting that the stream stars 
with different energies approximately follow a series of geomet- 
rically similar orbits, differing only in their time and length scales. 
The orbits of the stream particles predominantly sample radii in 
the range 5 < r < 150 kpc. Over this range, our best-fitting 
model potential for M3I is fairly close to a power law: $ oc 
with index k = —0.4 ± 0.2, significantly steeper than an isother- 
mal halo. Then we can use the concept of mechanical similarity 
(Landau & Lifschitz 1976) to relate different orbits that are similar 
in shape. Let the orbital time scale of such an orbit be proportional 
to a time scaling factor r, with t — 1 for the orbit of the progenitor. 
Then the radius scales as ^ (jjg velocity scales as -7-'=/(2-'=) ^ 

and the energy scales as T-2fe/(2-fe)^ while the kinetic energy T and 
total energy E relate as T = kE/{k + 2).To derive the phase 
space position of the stream, we consider a sequence of values of 
T. Let us choose the zero-point of time t to be the pericentric pas- 
sage where the stars are disrupted. (Our model explicitly assumes 
that the disruption is due to passage near M31, as opposed to some 
other cause such as an interaction with another satellite galaxy.) 
Given the orbital position fp{t) and velocity Vp{t) of the progeni- 
tor, we find the position of a stream star with orbital scale factor r 
asf = r^'''^~'''Vp(f/r), and the velocity as « = r*''*^"'''i7p(t/r). 
The tilt of the stream in phase space is then just a function of the 
time past pericenter. (If several streams exist from successive peri- 
centric passages, the temporal zero-point is different for each, im- 
plying a different scaling factor r for each at any given azimuth as 
well.) The extent of the stream, in contrast, depends on the spread 
in r or equivalently in orbital energy. This spread depends on the 
satellite mass as well as its orbit, so we postpone detailed consider- 
ation of the stream extent to the next section. With our model, the 
different orbits do not originate from precisely the same location, 
but the difference in starting location is smaller than the size of the 
progenitor itself so this is not a significant problem. 

To check this approximate model, we evolve an N-body model 
of the stream. We assume that the progenitor of the stream is a 
dwarf galaxy with total mass 2 x 10* Mq, consistent with the mass 
estimates from the dynamical properties of the stream of Font et all 
( 2005). We represent the satellite by a Plummer model, and take its 
scale radius to be 0.3 kpc (note the half-mass radius of a Plummer 
model is f» 1.3 scale radii). 

We use A'^ = 32768 particles in the satellite, with spline soft- 
ening length e = 30 pc; only about 0.1% of the mass is contained 
within one softening length of the center, indicating the satellite is 
well resolved. We initialize the particle distribution with the ZENO 
package of J. Barnes; the particle velocities are set by solving an 
Abel integral for the energy distribution (see lBinnev & Tremain3 
Il987ft . giving an initial configuration very close to equilibrium. We 
evolved this satellite in isolation for 2 Gyr, and found only negligi- 
ble changes in the structure of the satellite. 

We place the satellite on an inbound orbit on the near side 
of M3 1 . We take this to be the first approach of the satellite. In 
reality, unless the progenitor orbit consists of an extremely lucky 
shot, dynamical friction must be invoked to bring the satellite in to 
this point. 

We perform the calculation with the versatile parallel tree N- 
body code PKDGRAV (Stadel 2001; Wadslev, Stadel, & Oui^ 
2004). Our runs use the same fixed Hemquist bulge, flattened expo- 
nential disk, and Navarro-Frenk- White (NFW) halo potentials as in 
our orbit integration code above. The gravity tree uses an opening 
angle 9 — 0.8 and the node forces are expanded to hexadecapole 
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Figure 2. Test of the accuracy of our approximation for the position of the stream in phase space. In these panels, green dots show the particles in our N-body 
simulation of § 2.2, which uses our best fit to the southern stream data neglecting the stream-orbit tilt. The solid line shows the orbit calculated in this potential. 
Purple squares show the mean trajectory of the central particles of the satellite in this simulation, which closely follows the orbit. Red dashed fines show our 
approximation for the position of the stream. We choose the coordinate system for this figure so that the satellite orbit executes a radial loop within the X-Y 
plane, which thus also contains the stream. (Upper left): Spatial position of the particles, on an outbound portion of the orbit past the first pericentric passage. 
The axis ratio has been distorted for readability. (Upper right): Same as upper left, except shows the Y velocity on the horizontal axis. (Lower left): Same as 
upper left, except the satellite is now inbound. (Lower right): Same as lower left, except shows the Y velocity on the horizontal axis. 



order, while the individual particle timesteps are limited by the ac- 
celeration criterion Ai < ri{e/aY^'^ with 77 = 0.2. 

In Fig. |2| we demonstrate the accuracy of our approximation 
for the stream. To display the results more clearly, we choose a 
coordinate system just for this plot such that the radial loop after 
the first pericenter lies in the X-Y plane. The plots show the tra- 
jectory of a single test particle, and the median trajectory of the 
initially most-bound particles in the satellite. These nearly overlay 
each other, as expected; the small difference is due to the dynam- 
ical effect of the stripped outer portions on the inner core, as we 
can verify by decreasing the satellite mass. The upper and lower 
panels show the positions and velocities of the stream particles at 
two different time stages. At early times, there is a significant tilt 
of the stream in space relative to the orbit. Later on, this tilt is less- 
ee) 2005 RAS, MNRAS OOO.ITIfTTl 



ened, but the stream extends further in radius than the calculated 
orbit. This is a product of the systematic gradient in orbital energy 
along the stream, where the trailing particles have higher energy 
jjohnston. Sackett. & Bulloc3l200ll) . There is also a pronounced 
tilt in the space of Y versus vy- This relation is extremely tight, 
since the height reached depends almost entirely on the orbital en- 
ergy. Phrased another way, the large extent in physical space must 
be balanced by a small extent in velocity space in order to conserve 
phase space volume. 

Our approximation to the stream's phase-space location is 
shown as a dashed line. From here on we assume a potential slope 
of fc = —0.4; this is consistent with our power-law fit to the po- 
tential, and seems to produce the best results on this plot as well. 
Note the stream approximation describes only the central ridge of 
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the stream in phase space, and predicts neither its extent along this 
ridge nor its width. The stream approximation seems to be fairly ac- 
curate. Some deviations are apparent, but in general these are fairly 
small compared to the width of the stream. The approximation is 
certainly much better than assuming the stream follows the orbit. 
In the rest of the N-body simulations in this paper, we generally 
find good agreement with our stream approximation as well. The 
agreement worsens slightly as the satellite mass increases, since the 
self-gravity becomes more important, but it still appears acceptable 
given the observational errors. The approximation thus satisfies our 
need for an easy way to model the stream in a fitting routine, avoid- 
ing the need for many N-body simulations. 

2.4 Fits with stream-orbit tilt 

We now repeat the fitting exercise above, incorporating the stream- 
orbit tilt. This adds an extra parameter, the current phase of the pro- 
genitor along the orbit. We express this as a fraction of the radial 
period, i.e., Fp = t/tl^\ where t is the time since the progeni- 
tor went through the stream's last pericentric passage, and tr is 
the radial period of the orbit of the progenitor. (This is not to be 
confused with the azimuthal phase ^ employed by Johnston 1998 
among others, though the two are obviously related.) For a stream 
star, the orbital phase is instead F = t/tr, with tr the radial pe- 
riod of the particular star. The time scaling factor is then conve- 
niently given by r = tr/ti^^ = Fp/F. We keep y{0) fixed, so 
we now have eight free parameters: six orbital parameters, plus the 
halo parameter fh and the systematic shift Az in the depth. We test 
values for Fp ranging from 0.5 to 2.5, or in other words from the 
apocenter of the stream's radial loop to the one that is two radial 
periods ahead. We find Fp is nearly unconstrained by the fit; the 
changes it introduces in the stream location can be compensated 
for by changes in the other free parameters. We hence treat it as 
an external parameter, instead of trying to fit it. For the sake of 
discussion, we choose Fp = 1.0 in this section unless otherwise 
specified. 

When we allow fh to vary freely, without the inclusion of the 
galaxy fit term in , we find best fit values of /h > 4 or > 40. 
This is not only strongly disfavored by our galaxy fits but implies 
a halo concentration C200 ^ 7, which is at the bounds of physical 
plausibility (see Paper I). This fh value is higher than that in § 2.2, 
because the gradient of orbital energy along the stream adds to the 
required potential gradient. However, once we add the term from 
the galaxy fit, we find a best fit of fh ~ 1.2 ± 0.4 practically inde- 
pendent of Fp, similar to the results in § 2.2. This implies a mass 
M125 = (7.4± 1.2) X 10" Mq. The best fit orbital parameters for 
several choices of Fp are listed in Table|2|(orbits 2, 3 and 4). The 
reason for the lowering of fh from > 4 to ~ 1.2 is that the vari- 
ation in galactic dominates the variation from the stream. This 
is partly because the fit has another way to increase the computed 
potential gradient besides changing : Z\z is set to a lower value, 
placing the near end of the stream at a smaller distance from M3 1 . 
Recall that the systematic error in the depth is 16 kpc, compared to 
the error in the transverse position of about 2 kpc. We are roughly 
looking down the stream, so the distance uncertainty allows a large 
shift in the radius and thus in the potential. 

Even though the best-fitting fh is about the same as in § 2.2, 
the orbit is only about half as long as before; the radial period 
is 0.80 Gyr, apocenter is at 63 kpc, and pericenter is at 2.7 kpc. 
The stream-orbit tilt helps to explain why the stream is so strongly 
radial; the initial conditions required to produce large apocenter- 
pericenter ratios are rather less special than for an orbit that follows 



the stream. The apocenter and radial period decrease with Fp, since 
it controls the relative sizes of the stream and the progenitor orbit 
through the scaling factor r. For example, with Fp = 1.5, the best- 
fitting orbit has radial period 0.56 Gyr, apocenter 46 kpc, and peri- 
center 2.4 kpc. The period varies only slightly with fh - Typically, 
the projection of the second lobe on the sky lies at a smaller posi- 
tion angle, or northward of the NE side rather than eastward as in 
§ 2.2, but the direction is sensitive to exactly where the orbit passes 
through the disk. 

Using our best-fitting orbital parameters (with Fp — 1.0), we 
now repeat our N-body simulation. In this case we doubled the 
satellite mass to 4 x 10* Mq and raised the satellite scale radius to 
0.6 kpc, to get a better match to the length of the stream, and set the 
softening length to 0.1 of the new scale radius. We again start the 
satellite inbound towards the last pericentric passage of the stream; 
the initial conditions for this point are given as "2e" in Table |2| 
Fig.|3|compares our orbit and the resulting stream to the observa- 
tions. The agreement of the observations with the stream approxi- 
mation is excellent. The agreement with the N-body results is not 
quite as good, since as in Fig.|2|the approximation slightly underes- 
timates the azimuthal phase of the particles near pericenter, but this 
offset is comparable to the observational errors. As expected, the 
progenitor has a more compact orbit than the stream stars; the orbit 
turns around at about 2.5°, even though the stream itself extends 
out to 4.5°. With the chosen Fp = 1.0, the satellite lies against the 
M31 disk where it may have escaped detection. 

The stream's continuation on the opposite side of M31 is be- 
ginning to develop a shell feature in Fig.|3| This feature is not al- 
ways apparent; it depends on the orbit orientation and orbital phase. 
The continuation of the stream is also significantly rotated to the 
north compared to Fig. Q After the first pericentric passage, the 
tidally stripped debris gradually spreads out along the orbit, form- 
ing a well-defined stream. After the second passage, the stream ac- 
quires a significant width as well. This supports the idea that the 
narrow southern stream was indeed disrupted less than one radial 
period ago. The debris forward of the satellite is significantly dis- 
placed from the position of the orbit, just like the trailing debris. 
This effect should be taken into account when trying to match up 
possible stream orbits with the faint features in the M31 disk. 



3 STREAM AND PROGENITOR PROPERTIES 

In the previous section, we have considered the properties of the 
stream that have a direct bearing on the orbital parameters. We 
have obtained estimates of the orbit of the progenitor, but the or- 
bital phase of the progenitor is unconstrained by the fit and intro- 
duces an unwelcome degree of uncertainty to the orbit itself. We 
also found that the position of the next orbital loop could point in 
a variety of directions, depending sensitively on the orbital param- 
eters. So before we try to match the orbit to various observed fea- 
tures, let us consider what we can learn about the orbital phase and 
other properties of the progenitor from stream properties such as 
its length, width, velocity dispersion, and luminosity. These results 
should stand up regardless of the identification of the progenitor. 

We begin by generating an extended sample of N-body sim- 
ulations. We use three different orbits, namely 2e, 3e, and 4e in 
Table|2|which are the best fits for Fp = 1.0, 1.5, and 2.0 respec- 
tively. We again select a point on the inbound loop before peri- 
center as a starting point for each simulation. For each orbit, we 
test several values of the satellite mass. We adjust the scale length 
of the satellite along with its mass to keep the characteristic den- 
sity constant, ensuring similar amounts of satellite disruption in 
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Figure 3. Results accounting for a tilt between the stream and the orbit. Quantities shown in the panels are the same as in Fig.Q The orbit is shown with 
purple lines. The corresponding position of the stream calculated using our analytic approximation is shown with blue dashed hnes. We conduct an N-body 
simulation using initial conditions taken from this orbit; particles in this simulation are shown with red points. Observational results are shown with black 
symbols. The remnant of the progenitor is a small clump near the center of M31, in keeping with the orbital phase of Fp = 1.0. 



each run. The relation between the mass and velocity dispersion 
in these runs is consistent with that observed in dwarf galaxies 
(e.g., De Rijcke et al. 2005), suggesting this characteristic density 
is reasonable. In all other respects, the simulations are like those 
described before. The runs and some statistics derived from them 
are listed in Table|3| The run already displayed in Fig.|3|is "runl 1" 
in this table. 

We examine the resulting particle distributions at the timestep 
determined by Fp. One trend is obvious, from the positions of the 
particles in phase space: as the progenitor mass increases at a fixed 
orbit, the length, width, and velocity dispersion of the stream in- 
crease significantly. This suggests we can obtain an estimate of the 
progenitor mass by considering thes e propertie s of the stream. This 
work follows up the ideas in, .Font et alU2003) . but we improve the 
results by using our revised orbits and some modified formulae, and 
checking the results with our N-body simulations. 



3.1 Length and stellar mass of the stream 



The length of the stream is related directly to the spread of ener- 
gies of its stars. For a satel l ite of mass Ms and pericenter Rs in an 
isothermal halo, Ijohnstoni il998h found the energy range of tidal 
debris tends to be spread out over the range — 2e < AE < 2e, 
where e = Vt?(GMs/Vi^i?peri)^''^' Ms is the satellite mass, Rperi 
is the pericentric distance, and Vc is the circular velocity at pericen- 
ter. By examining our small sample of N-body simulations, we find 
results consistent with the mass scaling in this formula. However, 
we find a somewhat smaller spread: in our runs, 97% of the parti- 
cles have AE < 1.3e after the first close pass. (AE is distributed 
roughly symmetrically around 0, but we use a one-sided test here 
since in some simulation outputs the forward stream is affected 
by a second pass through peri center.) Aside from this change, the 
derivation o f lFontetalJiiooi) for the length of the stream may not 
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Table 3. Parameters of the N-body simulations used to estimate the stream quantities. Each run is labeled by an arbitrary integer The 
orbital phase of the progenitor is denoted by Fp, and the initial conditions corresponding to the orbit label are listed in Table l2l Ma is 
the total initial satellite mass, and is the scale length of the Plummer model. In each case we take the softening length to be O.laj,. 
The stream mass Matream is defined as the total mass of N-body particles within the phase space region outlined by Equations Islll II 
The fraction of stars in the southern stream still moving radially outward is given by f outbound ■ The values for the angular FWHM of the 
stream are derived from Gaussian fits with variable baselines; the value marked with a question mark differs significantly from the true 
width of the stream, owing to the mismatch between the narrow field and the large true width. Finally, avr gives the rms dispersion in the 
velocity relative to the observer of the stars in the sfi'eam region, computed by biiming the particles as in Fig. 171 and subtracting the mean 
velocity within the bin before finding the total dispersion. 
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be valid for large energy spreads or for non-isothermal halos. We 
hence derive the relation as follows. 

Our mechanical similarity argument in § 2.2 implies that the 
orbital energy E = Eq -\- AE depends on the time scaling factor 



r as E/Eq = 1 - AS/lSol = r 



2fc/(2-fc) 



-1/3 



(The energy 



Eq of the progenitor is negative, if we are using the approximate 
power-law form for the potential; hence r increases with AE). We 
then find 

|3 



Ma = 



Rperi \Eq\ I ^ -l/3\ 



(5) 



Here Tmax is the value of the time scaling factor at the end of the 
stream, and equals the ratio of the orbital phase of the progenitor 



2Fp. We obtain Rp^ri and Eq from the orbit 



that the stream extends just about to apocenter where the phase 

F = 0.5, so Tmax 

of the progenitor, and we take Vc = 250 km s~^ 

We can also derive a constr aint on the mass from the ot her end 
of the southern stream. Field 8 o f lMcConnachie et aljiioosh which 
is the closest field to M3 1 and the furthest along in the orbit. Since 
this field does not necessarily represent the end of the stream, we 
can only derive a minimum satellite mass, valid when the progenitor 
is assumed to lie within the visible extent of the stream: 



Ma > ^"^^ 



(6) 



We obtain tfs by matching the m coordinate of Field 8 to the cal- 
culated m coordinate of the stream as a function of r. 

We can now estimate the satellite mass required to produce the 
visible extent of the stream. We use the best fit parameters evalu- 
ated for a series of Fp, restricting the value of to the values 1.2 
(our best fit for most values of Fp) or 2.0 (~ 2cr higher). The pre- 
ferred mass from Equation|5|is shown as the black lines in Fig.|4] 
The estimates from the two values of fh do not differ too much, 
indicating our result for the required progenitor mass is fairly ro- 
bust. The main thing to observe is the strong correlation between 
the orbital phase and the mass. The primary reason for this is the 
dependence on Fp through Tmax- As the orbital phase increases for 
a. fixed orbit, the value of the mass must increase in order to create a 



stream out to apocenter. There are also changes in the orbital quan- 
tities Eo and Rperi as the best-fitting orbit changes as a function of 
Fp, but these are smaller effects and nearly cancel out in any case. 
The minimum satellite mass from the forward extension in Equa- 
tion |6| is also shown in the plot. This places only a weak limit on 
the mass and orbital phase, which will shortly be superceded. 

The estimate of the satellite mass from the stream length has 
significant uncertainty, since we are relying on an analytic scaling 
relation which has not been rigorously established. Also, the "end" 
of the stream is not very precisely defined; since the energy distri- 
bution of stream stars drops off rapidly but continously with AE, 
if the observationally apparent end of the stream extends to apoc- 
enter, there should be a small fraction of stars (~ 3%) that are still 
moving outwards towards apocenter. We can test the formula with 
our N-body simulations. Table|3|lists the fraction of the total stars 
moving outwards; this is clearly a very steep function of the pro- 
genitor mass for a given orbital phase. If we require this fraction 
to be 3% of the total stars in the progenitor, and interpolate the 
results in the table, we obtain an estimate of the progenitor mass. 
(This mass is not very sensitive to the adopted threshold, since the 
number of stars at large AE is falling rapidly with E). These esti- 
mates, shown as filled circles in Fig.|4] seem to be consistent with 
Equation |5| This is not too surprising, since we used these runs to 
calibrate the numerical constant in this equation, but it is reassuring 
nevertheless. 

The stellar luminosity in the stream gives us an alternate 
means of estimating the progenitor mass. The red dashed line in 
Fig.|4]shows the estimated stellar mass Mstmam. ^ 2.4 x 10* Mq 
of the stream (Ibataetal. 2001a). Allowing for dark matter, this 
gives a lower limit on the total mass. This immediately produces 
a tighter limit on the mass and orbital phase than the length of 
the stream alone. However, the fraction of stars that lie within the 
stream region can be quite small, and varies strongly with the mass 
and orbital phase. For the sake of argument, let us assume that 
the energy of the stream stars is distributed uniformly in the range 
— 1.3e < AE < l.Se. We can use this to derive the fraction of stars 
that lie within the allowed range tfs < t < Tmax- In turn, this 
raises the estimate of the mass of the progenitor to 
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1.5 

Orbital phase 

Figure 4. Shows estimates of the required mass Ms of the progenitor satel- 
lite as a function of orbital phase Fp . The solid and dotted rising black lines 
show the satelHte mass needed to make the stream extend to apocenter, de- 
rived from Eguationlsl These lines use the best fit orbital parameters for the 
two values of = 1.2 (solid) and 2.0 (dotted). The noise at high orbital 
phase corresponds to jiggles in the parameters from our fit. which are 
not very stable in this region. The black circles show the estimate of the 
mass from our suite of N-body simulations, obtained by interpolating the 
fraction of stars that ai'e still moving outwai'ds at the timestep correspond- 
ing to orbital phase Fp. The blue, rapidly declining lines at the left side 
of the figure show the minimum satellite mass needed to make the stream 
extend to its closest observed point (Field 8), again using the two cases 
ff^ = 1.2 (solid) and 2 .0 (dotted). The ho rizontal red dashed line shows the 
observed stellar mass hbataet alJbOOlaft . The red gently rising solid line 
shows the estimate of the progenitor mass needed to put this amount of stel- 
lar mass in the stream, as derived from EguationlTl The red squai'es show the 
con'esponding estimate of the mass from our suite of N-body simulations, 
obtained by interpolating the mass of particles found in the stream. The two 
data points at each orbital phase are slightly offset from each other for the 
purpose of clarity. Finally, the magenta dot-dashed line shows t he estimate 
of the progenitor's mass from the metallicity of the stream stars jFont et all 
l2005ft. 



-1/3 



Ms = 2Mstr 



-1/3 



.^-1/3 
' max 



(7) 



As long as the stream overlaps the observed length, out to Field 
8 (where the steeply rising black and falling blue curves cross), 
this equation gives more restrictive limits on the progenitor mass 
than the observed stellar mass alone. The results are shown as the 
solid red curve in Equation |4| This curve is a lower limit on the 
total satellite mass, since if the progenitor is partly composed of 
dark matter, or if it is not fully disrupted, an even larger mass is 
needed to match the observed stream luminosity. Hence, in Fig.|4] 
Ms and Fp should lie on the portion of the black curve that lies 
above the red curve. This implies Fp > 0.9 and Ms > 4 x 10** Mq. 
The square points in the figure are obtained by interpolating the 
progenitor mass AIs as a function of the mass Alstream of stars 
in the stream region; this is defined below in Equations Isl ll 1> and 
given in Table |3| The square points lie slightly above the curve, 
in part because the progenitor is not quite fully disrupted, and also 
because the stars are starting to leak out from the sides of our stream 



region as well as from the ends. This validates our argument that the 
red curve is a lower limit to the mass. 

The results in this plot are not extremely precise, but still sev- 
eral points are evident. First, the progenitor most likely lies ahead 
of the stream in its orbit. Second, the required mass is strongly cor- 
related with the orbital phase. Third, the lower limit on the mass 
from the stream luminosity roughly matches the estimate of the 
mass from the stream length. This suggests that there was not much 
dark matter before last pericenter, and the progenitor was largely 
disrupted at last pericenter. It is quite plausible that the outer dark- 
matter-dominated regions were stripped by tides before the stellar 
core was disrupted. Fourth, the coincidence of these two curves is 
only obtained if we assume the stream is young (stripped at last 
pericenter). If we assume instead that the stream resulted from two 
pericentric passages ago, the estimated mass from the stream length 
(black curve) shifts down by a large factor, but the lower limit on 
the mass from the stream's luminosity (red curve) does not change 
much. In other words, the total mass implied by the directly ob- 
served mass in stream stars would be enough to spread out such an 
old stream to an unacceptable degree. Unless our estimate for the 
orbit is badly in error, we can rule out this model. 

Using the stream's metallicitv. lFont et alj J200^ estimated the 
progenitor had a stellar mass of Ms = 5 x 10^ ^Q- This estimate 
has an uncertainty of perhaps a factor 3, given the scatter in the 
metallicity-luminosity relation. Thus it is consistent with our sim- 
ilarly rough estimates of the satellite mass for almost any orbital 
phase, but the agreement is best for Fp '-^ 2.5. Alternatively, the 
progenitor may have lost stars to gravitational tides earlier in its 
life, moving its mass Ms at last pericenter off the mass-metallicity 
relation. In this case we expect Fp < 2.5. 

3.2 Other stream measures 

We now consider other observational measures of the stream com- 
pared to our N-body simulations. Our set of N-body models shows 
particles at a wide range of transverse distances from the center 
of the stream, not all of which would be picked up in the current 
observational surveys. The extent of this satellite debris depends on 
the mass and radius of the satellite as well as the pericenter distance 
Rp. To exclude the more extreme stellar trajectories in the stream, 
we define a "stream star" as one lying in the region 



0.60° 
-0.14° 
-30kpc 
-410kms~^ 



< 
< 
< 
< 



m < 4.70° 
n < 0.62° 
d' - 20m - 



Az <60 kpc 



70m < -210 km s 



(8) 
(9) 
(10) 
(11) 



The length and transverse width of this boundary approximately 
match the fields in McConnachie et al. (2003). The Az term com- 
pensates f or the fitted offset between our distance calibration and 
that of .McConnachie et al. (.2002). We chose a more generous se- 
lect ion boundary in velocity sp ace than used bv llbataetal1<2004 
and lGuhathakurta et alj i200 3). because we do not have the prob- 
lem of excluding a background density of stars; in any event the 
velocity-position relationship is tight enough that this is not a ma- 
jor factor. 

Fig. |5| shows the distribution of N-body particles along the 
stream. This is analogous to fig. 4 of^ McConnachie et alJi2003l) . In 
comparing the two, we should remember that "background" (non- 
stream) stars probably contribute a substantial and spatially varying 
baseline to the observed histogram. By monitoring the behavior of 
the histogram for different timesteps and different runs, we see that 
this plot is essentially another way to check the length of the stream. 
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Figure 5. Shows the distribution of N-body particles along the stream, 
using run 24. Position m along the stream is define d the same way as in 
Fig.|3](so its sign is reversed from Fig. 4 of McConn achie et alj2003ft . Cuts 
are made to the sample as described in the text. The dotted line continues 
the histogram without the sample cut in Eguationlsl 
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Figure 6. Shows the distribution of N-body particles as a function of n, the 
coordinate direction perpendicular to the stream, using run 24. Position n 
across the stream is defined the same way as in Fig.|3](so the sign is reversed 
compared to fig. 3 of McConnachie et al. 2003). The NE side of the stream 
is to the left. Cuts are made to the sample as described in the text. The low- 
amphtude magenta, red, orange, and yellow histograms (highest to lowest at 
the peak) show the distribution in quadrants of increasing distance along the 
stream, while the black histogram shows the total. The blue dashed curve 
shows a Gaussian fit with a variable baseline, having FWHM 0.57°. 



If the stream is too long, the histogram will rise towards its end; 
whereas if the stream is too short, it will cut off in the middle of its 
range. To the limited extent we can tell, given the background issue, 
the distributions appear to agree quite well. This again indicates our 
progenitor mass and pericenter are reasonable. 

In Fig.|S| we display the distribution of N-body particles in the 
direction normal to the stream. This can be compared to fig. 3 of 
pvlcConnachie et al. (2003), though again that figure contains back- 
ground stars not present in our run. The stream in our run happens 
to be centered in our chosen selection region (the observed stream is 
slightly offset from the selection region used here). We find the rise 
in the number on the NE side is comparable to that in the observed 
stream (the other side is cut off by the field boundary in the observa- 
tions). However, some outlier particles are present in non-Gaussian 
tails outside the width of the fields, as can be seen in the top panels 
of Fig.|3| Stars showing this behavior would probably fade into the 
general halo population in surveys using number counts alone. In 
fact, it is possible that they make some contribution to the wide- 
field, minor-a xis photometric surveys o dPurrell. Harris. & Pritchel 
<200lLl2004) . These fields are far enough away from the stream 
that, in our "run24", ^ 5 x 10^ Mq of stream material is covered 
by the two fields in this survey. From the relative surface bright- 
nesses in the stream and in the survey fields, we estimate this rep- 
resents only a few percent of the total mass of halo stars in these 
fields. However, this contamination would be very sensitive to the 
orbit, size, and radial mass profile of the incoming progenitor, so we 
cannot rule out the possibility that it contributes significantly to the 
observed population in these surveys. Interestingly, these surveys 
found a metallicity distribution peaking at [m/H] ~ —0.5, similar 
to the metallicity of the stream. In general, these outlier stars will 
share in the distance and velocity signature of the stream, so they 
may be reco gnizable in spectroscopic surveys. 

iFont et al. (2005.) suggested that the spatial width of the stream 
(defined to contain 80% of the stream luminosity, which for a Gaus- 
sian is the FWHM) at a physical distance R from M31 could 
be expressed as lu = sR, where s = {GMs /V^ RperiY^^ ■ 
Using this estimate, they fitted a Gaussian with a floating base- 
line to the transverse distribution of the stream stars and inferred 
Ms ~ 1.3 X 10* M!©; this is quite low compared to their estimate 
from the stellar metallicity. Font et al. dismissed this discrepancy, 
noting that if the progenitor was partially supported by rotation it 
might reduce the width of the stream. 

We mimic their procedure by fitting a Gaussian with an ar- 
bitrary baseline to the N-body runs. The resulting estimate for the 
FWHM is given in Table|3| which shows that it increases system- 
atically with Ms. We find that this is not a completely robust pro- 
cedure to recover the true width of the stream, because the window 
is narrow and the non-Gaussian tails can masquerade as a baseline. 
Nevertheless, in most cases we recover nearly the same answer as 
we would with twice as wide a window, so their width measure- 
ment is probably accurate. Font et al.'s analytic formula systemati- 
cally overpredicts the estimated width in our runs. However, if we 
multiply their formula for w by 0.4, we get fairly good agreement 
overall. It may be that their formula simply needs this scaling con- 
stant, which raises the mass at a given width by more than an order 
of magnitude, although more thorough testing of this is warranted. 
The best agreement of our N-body runs with the observed width is 
obtained for Ms ^ 10^ consistent with our other estimates. 
The mass estimated this way is not very precise but is probably 
good to a factor 3. If we interpolate the results in Table |3| to find 
the FWHM at a stream mass of Mstream — 2 x 10** Mq, as we 
did with the mass earlier, we find results of 0.71°, 0.50°, and 0.70° 
for Fp — 1.0, 1.5, and 2.0 respectively. Although the precision of 



© 2005 RAS, MNRAS 000,[T][ni 



Investigating the Andromeda Stream -II 13 



12 3 4 

Distance m along stream (degrees) 



runs occupy only a small range in azimuth, blurring any correlation 
that might exist. 

If we again interpolate the results in Table |3| this time to 
find the velocity dispersion at a stream mass of Mstream = 2 x 
10** Mq, we find results of 18, 18, and 20kms"^ for Fp = 1.0, 
1.5, and 2.0 respectively. Again, there does not seem to be a strong 
trend with orbital phase for a. fixed mass of stars in the stream. 

In summary, we find that the various observed properties of 
the stream give consistent indications of the mass. Our simulation 
with a mass of Ms — lO'"* Mq and an orbital phase of Fp — 1.5 
is the one that is most consistent with the observed properties of 
the stream. There is some leeway in these quantities, however, as 
long as the mass and orbital phase increase together in keeping with 
Fig- 111 Aside from detection of the progenitor itself, radial veloc- 
ity surveys will probably be the most powerful tool to constrain 
the progenitor mass. Radial velocity surveys to a greater transverse 
extent may help to differentiate between unrelated halo stars and 
the sideways extension of the stream, which would give additional 
constraints on the properties of the progenitor. 



Figure 7. Shows the line-of-sight velocity dispersion of N-body particles 
in run 24 (red squares), as a function of position m along the stream. The 
points are computed by taking particles within the stream region (Equa- 
tions |8]^3' binning them in m, and finding the standard deviation of the 
velocity of the particles in each bin. Black circl es with error bars show the 
observed velocity dispersion in five fields from lGuhathakurta et alJ 
andfperauson et al. ( 2004). We indicate the field in Ferguson et al. closest 
to M31 with an empty symbol due to its potential contamination with the 
disk. A rough error bar of 10 km s~^ is assumed for each observed point. 



these estimates is low, they indicate that models that reproduce the 
stream mass also reproduce the stream width fairly well, indepen- 
dent of the orbital phase. 

Fig. Q shows the rms velocity dispersion at several points 
along the stream. Observed values of the velocity dispersion are 
shown in the figure as well; these are quite uncertain, given the 
small number of stars and the possibility of contamination by ad- 
ditional structure in the background density. The dispersion in the 
N-body model is certainly consistent with the observations, sug- 
gest ing our progenitor mass is roughly correct. 

iFontet aljT2005h suggested > 2.0 x 10* Mq from the ve- 
locity dispersion; they noted that this is a lower limit since the ob- 
served portions of the stream may be dynamically colder than the 
progenitor itself, which we indeed find to be the case (cf. the right- 
hand panels of Fig. |2}- We find the central velocity dispersion of 
the progenitor is ~ 25 kms~^, while the M31 -centric radial veloc- 
ity dispersion as a function of M31 -centric radius drops to very low 
values ( ^ 5 km s~^) throughout the stream. The larger line-of-sight 
dispersion in Fig.Qis primarily due to the fact that the stream stars 
form a fattened, tilted tube in three dimensions, and thus a single 
line of sight probes particles at a range in M31 -centric radius which 
depends on the tube's thickness. As a result, there is a good corre- 
lation between the transverse width and the velocity dispersion in 
Tabled 

Font et al. also suggested that the velocity dispersion, should 
increase sharply near apocenter. In contrast. Fig. shows it to be 
roughly constant to the last measured point, at nearly the end of the 
stream. The difference may arise from Font et al.'s use of azimuthal 
binning rather than our binning in transverse position m. It is dif- 
ficult to compare our results directly, as the stream particles in our 



4 FORWARD CONTINUATION OF THE STREAM 

From the analysis in the previous section, we now have a better 
idea of the mass and orbital phase of the progenitor. We now con- 
sider the continuation of the stream past the next pericenter. For 
the fit shown in Fig.|3| the second lobe projects about 30° counter- 
clockwise of MSl's NE major axis, but this direction is sensitive 
to the orbital parameters. If we impose further observational con- 
straints, this loop can point in a variety of directions which encom- 
pass a number of faint features observed near M31. As photometric 
and velocity surveys go deeper and deeper, more faint features will 
probably continue to turn up, and it is quite possible that the pro- 
genitor is not even a currently identified feature. Nevertheless, it is 
worth testing which of the currently identified features can be the 
progenitor. 

Table|4|gives the observational data for a number of these fea- 
tures in the disk and inner halo regions of M31. We perform fits 
in the same manner as in § 2.4, but we add further data points for 
the sky position of the progenitor, and its depth or radial veloc- 
ity if these are known. Most of these possible progenitors require 
a phase Fp > 1.0. Once the progenitor passes through pericenter 
for a second time, a second stream develops, which again does not 
follow the orbit. However, in our N-body simulations, the progeni- 
tor itself follows the calculated orbit fairly closely, though there are 
some slight deviations due to the self-gravity of the satellite matter. 
The position and velocity of the stream thus can be derived directly 
from the calculated orbit and the orbital phase Fp. 

The errors on the observational data are in many cases not well 
specified, so we have chosen reasonable values in order to perform 
the fit. For example, some of the features are not very well-defined 
spatially, including the Northern Spur, the Merrett et al. planetary 
nebulae (PNe), and the Eastern Shelf. We assume positional er- 
rors of 0.05-0.2°, depending on how well-defined the object is. 
In Table|4| if the distance or velocity are not known for a particu- 
lar feature, we supply the results from our fit in brackets. In some 
cases there may be alternate solutions to the one we have com- 
puted, where the progenitor has a different orbital phase, especially 
if only the sky coordinate of the feature is known. We now discuss 
the suitability of the follow ing features as p rogenitors. 

Linear conft'«Mario«: iMcConnachie et aL 12003.) observed two 
fields along the linear projection of the southern stream on the other 
side of the disk. Even though there is no obvious projection seen in 
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Table 4. Kinematic data and predictions for possible progenitors of the stream. Brackets indicate tliat the field is a prediction from 
our best-fitting orbit, as opposed to observational data, d' and give the radial displacement and velocity relative to M3 1 . The depth 
offset Az and orbital phase Fp are taken from our fit. The final entry gives the value of reduced from our fit. The number of degrees 
of freedom in each case is 13 plus the number of observational constraints in the table. 
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Figure 8. Best-fitting orbit passing through M32, and the particle distribu- 
tion in the coiresponding N-body simulation. Crosses, orbits, and dots have 
the same meaning as in Fig. The particles con'esponding to M32 itself 
are barely visible in this plot as the clump just south of the center of M31. 
The point density in the stream is smaller than in Fig. |3] because we use 
a larger particle mass in this run to match the larger mass of the progeni- 
tor. Observed positions of some of the faint features are plotted as symbols: 
Northern Spur (cross), And NE (circle), eastern Merrett PNe (triangle), the 
Eastern Shelf (square), M32 (star). And VIII (diamond), and the Gl clump 
(star). 



surface brightness maps, they obtained red giant color-magnitude 
diagrams for these fields consistent with a stream component. We 
assign a position to this feature using the center of Field 13 and 
errors of 0.2°. Since the projected orbit tends to bend substantially 
at pericenter, we find it is very difficult to fit this feature as the 
prog enitor itself, in agreement with the conclusions of I bata et alj 
i2004i,') an d-Font et aL 1,2005.) . However, it is possible that the stream 
scatters enough debris to the side when it passes through pericenter 
to produce the stars seen by IVIcConnachie et al. 

NGC 205 and nearby arc: NGC 205 is close to the projection 
of the stream to the other side of ]VI3 1 . However, its velocity has 
the opposite sign from that of the stream. If it is the progenitor, 
it must lie nearly at the next pericenter, i.e. at phase ~ 2.0. 
Because of its location north of M31, it is difficult to fit an orbit 




Xm3i (degrees) 

Figure 9. Best-fitting orbit passing through the Merrett et al. (2003, 2004) 
PNe. In this plot we have changed to their M31 -aligned coordinate system 
to allow easier comparison with their figures. X^'n increases along M31's 
major axis to the SW, and Ymh along its minor axis to the NE. The four 
crosses in each panel show approxi mate points along the narrow-velocity- 
dispersion group of PNe identified in lMerrett et alj i2004ft . 
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with this constraint. The st ellar arc through the center of NGC 205 
iMcConnachie et all2004l h produces the same difficulties with the 
orbit as NGC 205 itself. In addition, the bluer color of the red giants 
in NGC 205 and the arc indicate a lower metallicity compared to 
the stream, again arguing against a connection with the stream. 

M32: M3rs other prominent nearby satellite galaxy is an ob- 
vious candidate for the origin of the stream, and was proposed as 
such as soon as the stream was discovered ( Ibata et al. 200 Ig)- M32 
lies almost directly on the path of the stream. However, as with 
NGC 205, its radial velocity has the wrong sign for it to lie within 
the southern stream. Hence it too would have to be at the next peri- 
center (Fp ~ 2.0). Interestingly, on orbital grounds alone, it is 
not possible to rule out this possibility. We find an acceptable fit, 
even though the observational constraints on this fit are more nu- 
merous and more precise than for any other candidate (since the 
position, distance, and radial velocity are all fairly well known). 
The angle of the second loop in this fit is not very robust, and 
could accomodate a wide range of position angles ranging from 
to 60°, which encompasses a range of observed structures, al- 
though the best-fitting orbit, shown in Fig. |8| does not result in 
an obvious match to any in particular. The progenitor mass before 
last pericenter required to produce the length of the stream is about 
Ms — 2.5 X 10^ Mq. However, it then must lose much of its 
mass to the forward and backward stream s in or der to produce the 
observed luminosity of the stream. iMateol il998^ gives the current 
mass of M32 as 2.1 x 10^ Mq. Within the uncertainties, the masses 
are p robably co nsistent. From studying the isophote ellipticity of 
M32. IChoi.'Guhathakurta, & lo hnston ( 20QZ) inferred that it is on 
a highly eccentric orbit with JJperi < 1.7kpc, which is consistent 
with our best-fitting orbit. However, these authors also inferred that 
M32 is most likely moving eastward; this is not consistent with our 
orbit, though it is difficult to evaluate the reliability of their argu- 
ment. There are other difficulties with our specific orbit: it creates 
a plume in the northward direction that has not been reported ob- 
servationally, and it is not clear that the stream forward of M32 can 
match the observed structure. The metallicity distribution of M32 
is complex, and it is not clear whether or not it would be consistent 
with that of the stream. Nevertheless, on balance of the evidence, 
M32 continues to be a p ossibility for the prog enitor of the stream. 

Merrett et al. PA^eil Merrett et al] i2003h surveyed the disk of 
M3 1 for PNe and found a number that had velocities inconsistent 
with the disk kinematics. Almost by definition, most of them move 
in a direction opposite to the local rotation of the disk, though a few 
move in the same direction but with larger speeds. Merrett et al. 
i2004l pointed out a subset of these on the northeast side that oc- 
cupy a narrow region in position- velocity space. We can easily fit 
this group as the progenitor of the stream. (Merrett et al. also fit 
an orbit connecting this group to the stream using their simpler 
potential, but without demonstrating agreement with the detailed 
properties of the stream.) Choosing the position and velocity to lie 
within the northeast group, we find an excellent fit with the initial 
conditions given in Table|2| Fig.|9|shows the positions and veloc- 
ities obtained from an N-body simulation using this orbit, using a 
coordinate system aligned along the "^^^j"'" of 1 • These can 
be compared to fig. 3 of Merrett et al. '2004V Interestingly, there 
are a significant number of N-body particles occupying larger ve- 
locities at a given position than the main concentrat ion, but none 
with s maller velocities; a similar pattern is seen in iMerrett et alJ 
i2004) . This feature is thus a prime candidate for the progenitor of 
the stream. 

N orthern Spur. The star-count maps of M3 1 in lFerguson et alJ 
i2002h show a faint clump of stars at the NE end of M31, called 
the Northern Spur. We estimate the position of the Spur from the 
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Figure 10. Best-fitting orbit passing through the Northern Spur, and the 
particle distribution in the corresponding N-body simulation. Crosses, or- 
bits, and dots have the same meaning as in Fig.|3] Object symbols are the 
same as in Fig.|8] The progenitor is visible as the dense elongated clump 
near the NE end of M3 1 . 



image in I Zucker et alJ i2004h . which also is suggestive of a loop 
of emission extending from the Spur. The red giant branc h color 
is similar to that of that of the stream (Ferguson et al. 2002). This 
feature is a promine nt candidate for the continuation of the stream 
jMerrett et all2003h . 

Some radial velocity observations exist in the Spur region. 
From spectroscopy of red giants, Ferguson et al. 1 2004) claimed a 
velocity of ~ 150kms~^, significantly different fro m that of the 
disk in this region, suggesting it is a distinct object. iMerrett et all 
( 2004) found an overdensity of PNe in this region. The velocities 
of some of these are consistent with the disk velocity in this region, 
some have larger velocities, and some have the wrong sign alto- 
gether. These authors inferred that the Spur represents a warp in 
the disk of M3 1 . The possible contamination from the disk and the 
presence of velocities of both signs suggest that if the progenitor 
lies in the Spur region, its velocity cannot be reliably constrained at 
present. 

If we impose no velocity constraints, we find an excellent fit to 
the position of this feature, as seen in Fig. llOl The predicted loca- 
tion is on the near side of M3 1 , and the progenitor is outward bound 
from pericenter, giving a negative radial velocity. Thus slight dif- 
ferences in its color-magnitude diagram in Ferguson et al . (20oJ), 
as compared to that of the stream, could possibly be explained by 
its smaller distance. One possible problem is the faint luminos- 
ity of the Spur ; from the area and surface brightess discussed in 
iFerguson et alj j2002), we estimate its V-band luminosity is only 
8 X 10® Lq, less than that of the stream. The simulation shown in 
the figure produces a fairly dense group of particles at the position 
of the Spur, and would probably need to be adapted to produce a 
fainter, lower-surface-brightness core. 

If we instead assume that the velocity is positive, as suggested 
by the majority of the velocity data, the progenitor must be moving 
inwards once again. In this case it is more difficult to fit a model to 
the Spur. Our best fit is formally ruled out, though this conclusion 



© 2005 RAS, MNRAS OQO.ITIUtI 



16 M. A. Fardal et al. 



is dependent on the assumed errors for the Spur's properties (we 
assume a velocity 200 ± 30 km s~^ and a positional error of 0.2°). 

And NE: This low surface brightness enhancement far ou t 
along the major axis of M31 was discussed by Zucker et al. '('2004"). 
We find it is quite easy to make the orbit pass through this object. 
However, this solution is objectionable on other grounds: there is 
no visible connection to the stream, despite the low surface bright- 
ness of M3 1 at these projected radii, and the luminosity is probably 
too low for it to be an intact progenitor. 

Eastern shelf: Iperguson et alj J2002l) discussed a diffuse sur- 
face brightness enhancement on the eastern side of M31, which 
contained high-metallicity red giants judging by their color. This 
feature is now known as the "Eastern Shelf". This color may be 
consistent with that of the Northern Spur and of the stream itself. 
With the position we have assigned to this feature, it is not easy to 
find an orbit where this feature is the progenitor itself, despite the 
weakness of the observational constraints. This conclusion differs 
from that in Font et al. (2005) because the greater strength of their 
potential, and their use of an orbit that follows the stream, affect 
the direction of the second loop. The low surface brightness of this 
feature also argues against it being the progenitor. It might repre- 
sent debris from the continuation of the stream, though for none of 
the orbits we tried did the debris reach the position angle of this 
feature. 

And VIII: This is another group of objects with abnormal ra- 
dial velocities in the M31 disk (Morrison et al. 2003), in this case 
a combination of PNe, globular clusters, and H I detections. This 
feature lies nearly across the stream and has a velocity consistent 
with it, but it is extended for ~ 1° transverse to the stream, which is 
not expected for a progenitor moving along the stream. Hence if it 
is the progenitor, it has probably passed through pericenter already. 
Incorporating this condition as a constraint, we find that it is very 
difficult to fit an orbit to this fea t ure. 

GI Clump: iFerguson et alj i2002h identified a diffuse surface 
brightness enhancement in the vicinity of the globular cluster GI, 
which they called the "GI clump". We did not attempt to fit this as 
the progenitor of the southern stream. By virtue of its position, this 
would have to lie at i^p ~ 2.5. This would make for a very mas- 
sive progenitor. The GI clump, however, is less luminous and thus 
less massive than the stream, and it is not clear that the progenitor 
mass in excess of these two structures can be hidden. In addition, 
the color-magnitude diagram of the clump differs from that of the 
southern stream. However, it is at least suggestive that in a number 
of our N-body simulations of other possible progenitors, there is a 
clump or shell-like feature in roughly the vicinity of this feature, 
where the forward section of the stream turns around (see Figs.|8| 
andfTol. 

In summary, the preferred mass and orbital phase of the pro- 
geiutor indicate there is probably a significant amount of structure 
near M3 1 that is physically connected to the stream but ahead of 
it in its orbit, amounting to > 10® Mq in total. We can probably 
rule out most of the objects on the southwestern half of M31 as 
progenitors, with the possible exception of M32. We can fit rea- 
sonable orbits through several features in the northeastern half: the 
Northern Spur, And NE, and the Merrett et al. PNe. N-body simu- 
lations following these orbits produce a spray of debris as well as 
a number of discrete-looking features; there is some resemblance 
between this debris and the structure seen in the observations, even 
though we do not have a detailed match. We have made significant 
approximations in the fitting and simulation process, among them 
the neglect of dynamical friction, the assumption of a static fixed 



potential, and the neglect of the satellite's self-gravity in deriving 
the best-fitting orbit. Thus we have not tried to refine the match to 
the observed features any further. However, by predicting the dis- 
tance and velocity of each feature, we hope to speed the process 
of confirming or rejecting a connection between the various faint 
features and the southern stream. 



5 SUMMARY 

In this paper and Paper I, we have attempted to construct a model 
for the interaction of M3 1 and the southern stream. We constructed 
a simple analytic model for the potential, and found its best fit pa- 
rameters by fitting to a number of observational constraints. We 
then combined this galaxy fit with a fit to the observed path of the 
southern stream. We found the best fit parameters as a function of 
the orbital phase of the progenitor. By considering the effect of the 
progenitor's mass on the observed properties of the stream, we de- 
rived an N-body model of the interaction that fits current observa- 
tions reasonably well. We tested the orbital fits against various ob- 
jects and surface brightness features near M31, ruling out some as 
the progenitor and making predictions for the kinematic properties 
of the others. 

In comparison to the stream from the Sagittarius dwarf galaxy 
that surrounds the Milky Way, there is currently much more lee- 
way in fitting the Andromeda stream. The progenitor of the former 
stream is known, and its proper motion has been measured. Dis- 
tance errors for this stream are a factor of 10 smaller than for M31. 
In addition, this stream has been followed for nearly a full circle 
around the Galaxy, as opposed to the Andromeda stream which 
spans less than half a radial period. Thus, this stream is currently 
a much better probe of its parent halo's properties than in the case 
of M3 1 and its southern stream. In this paper, we have found that 
the orbital phase of the progenitor is a major source of uncertainty 
in the orbital parameters. If we can identify the progenitor location, 
we can immediately put much stronger constraints on the orbital 
properties of the stream. We can then get a better understanding of 
the mass distribution in M3I, calculate the properties of the pro- 
genitor satellite, untangle the possible relationships of the various 
surface brightness features, and perhaps find signatures of the pass- 
ing satellite in the M3 1 disk. 

Our specific conclusions are as follows: 

• The stream follows a path in phase space significantly differ- 
ent from the orbit of the progenitor satellite. Allowing for this ef- 
fect decreases the estimated apocenter and increases the estimated 
gravitational force in the M31 halo. 

• The narrowness and surface brightness are consistent with a 
young stream, resulting from the last pericentric passage of its pro- 
genitor. Larger ages are ruled out by the narrowness and high lumi- 
nosity of the stream. 

• With the current observational uncertainties, the stream itself 
does not add much constraining power to the galaxy model. Our 
best combined fit for the mass of M31 within 125 kpc is (7.4 ± 
1.2) X 10^^ Mq, only slightly larger than in our best-fitting model 
from Paper I. 

• If the progenitor satellite has orbital phase Fp — 1.0 (at the 
next pericenter beyond the stream location), our best fit to the orbit 
of the progenitor has radial period 0.8 Gyr, apocenter 63 kpc, and 
pericenter 2.8 kpc. However, Fp, which determines the size and 
time scales of the orbit, is poorly determined. We find it likely that 
the satellite or its disrupted remnants lie beyond next pericenter. 

• The dynamical mass of the progenitor at the previous pericen- 
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trie passage was Ma « 10^ Mq, within a factor 3. Its mass was not 
heavily dominated by dark matter; otherwise, the surface brightness 
of the stream would be less than is observed. 

• Depending on the exact properties of the progenitor, signifi- 
cant amounts of debris from the progenitor may lie outside the cur- 
rently observed width of the stream. These stars could contaminate 
the metallicity distribution measured in off-stream fields, but can 
be detected since they share the velocity and distance signatures of 
the stream. 

• Most of our orbital fits would place the continuation of the 
stream past next pericenter and to the NE of MSl's center. We can 
fit reasonable orbits through several observed features, among them 
the Northern Spur, And NE, M31, and a feature in the planetary 
nebula distribution found by Merrett et al. The latter feature seems 
to be the most plausible candidate for the progenitor of the stream. 
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